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Consider "lagged" Fibonacci sequences a(n) = a(n — 1) + a([n/k\) for k > 1. We show that 
lim n _> 00 a(kn)/a(n) • In n/n = kink and we demonstrate the slow numerical convergence to this 
limit and how to deal with this slow convergence. We also discuss the connection between two 
classical results of N.G. de Bruijn and K. Mahler on the asymptotics of a(n). 



I. INTRODUCTION 

Let k > 1 be an integer and consider "lagged" Fibonacci type sequences 

Ofe(n) = a k (n - 1) + a k ( ^ ) (1) 

with initial value 

a k (0) = l. (2) 

These "almost linear recurrence" has many interesting arithmetical properties [1]. The value a k (n) equals the number 
of fc-ary partitions of kn, and the corresponding sequences are listed in the OEIS as A000123 , A005704 , A005705 and 
A005706 for k = 2,3, 4, 5, respectively. In this contribution we will study the asymptotical behavior of the ratio 

a k {kn) Inn 

Ck(n) = . (3) 

a k (n) n 

The OEIS entry for A000123 quotes a conjecture due to Benoit Cloitre, claiming that 

lim Cfe(n) = const. = 1.63 .... (4) 

n — >oo 

The same conjecture (but with const. = 1.64 .. .) appears for the related sequence A033485 . We will prove that the 
essential part of the conjecture (existence of the limit) is true, but that its numerical part is incorrect. In particular, 
we will apply a classical result of de Bruijn [2] to prove that 

lim Cfe(n) = kink. (5) 

n — >oo 

Note that 2 In 2 = 1.386 . . ., which differs significantly from the value in (4). 

In the second part we will discuss the rate of convergence of c k (n). It turns out that this rate is so slow that 
straightforward numerical measurements of c k (n) cannot be used for an accurate measurement of c k (co). This may 
explain the inaccurate numerical value in (4). It turns out that another classical result on the asymptotics of a k (n) 
due to K. Mahler [3] can be used as a device for an accurate numerical determination of c k (n) all the way to the 
asymptotic regime. 

In the final part we will discuss the connection between the two asymptotic formulas of de Bruijn and Mahler. 
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II. ASYMPTOTICS 



Using an integral representation (Mellin transformation) of the generating function for afc(n) and a saddle point 
integration, dc Bruijn [2] showed that 
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lna fcW =^rr( ln » - In Inn) + - + — + — — Inn - 1 
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where ipk is a periodic function with period 1, 



j = -oo 



The Fourier coefficients are 



and 



1 r f^ij\ 



2nij 
Ink 



(J * 0) 



a (fc) 
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-7i - ^7 2 + — 7r 2 + — In 2 fc 
' 2 ' 12 12 



(7) 



(8) 



(9) 



where 7 = 0.5772156649 ... is the Euler constant and 71 = —0.0728158454 ... is the first Stieltjes constant. 

The Fourier series for ?/>fc(x) converges absolutely and uniformely because the coefficients ocj(k) decay fast enough: 



and 



|r(«)| = o(\t\-h~^) 

\<;(l+it)\=OQn\t\). 



Plugging (6) into (3) provides us with 

lnc fc (n) = \n(k\nk) + A^fc(n) + O 

where 



(In Inn) 
Inn 



AMn) = *> ( log * (icih) - logfe ( x + 1^)) " *" ( logfc (1^)) 



(10) 

(ii) 

(12) 
(13) 



Intuitively Atpk(n) should vanish for n — > 00, but to be sure we need to investigate the Fourier series for ipk in more 
detail. In particular, we have 



|A^ fc (n)| 



J2 cxp (2wij log fc (j^y ) (exp ^-2«jlog fe (l + ^ 



j = -oo 



- 1 



cxp -2m j log fc 1 



1 



1 



logfc n 



logfc 1 



!ogfcn 



j=-oo 

CX) 

In the last line we have used the inequality 

|e ix -l|<M (ieR). 
Now because of (10) and (11) we know that ■ \jaj(k)\ < 00, and hence 



A^(n) =0(1/ logfc n). 



(14) 
(15) 



This concludes our proof of (5). 




Figure 1: Numerical evaluation of C2(n) together with &2(fi) from (19). The difference between C2 and hi is smaller then the 
linewidth for all n > 100. The dotted vertical lines indicate n = 10 3 , 10 6 and 10 9 . 



III. NUMERICAL EVALUATION 



The recurrence (1) appears in the analysis of the Karmarkar-Karp differencing algorithm for number partitioning 
[4] . In this context we learned that the convergence to the asymptotic regime can be extremely slow. We will see that 
this is also true when we try to probe the asymptotics of c k {n) numerically. 

To calculate c k (n), we need a(kn) and a(n), but because of 

n-l 

a k {kn) = Ofc(n) + k ^a k (j) , (16) 
j=o 

the value of a k (n) plus the sum of the preceeding terms is sufficient. The bottleneck for calculating a k (n) is memory, 
not CPU time, since n(l — 1/k) values must be stored to compute ak(n). We used the Chinese Remainder Theorem 
to keep the individual numbers small and managed to calculate C2(n) for n up to 3 • 10 9 on a PC with 4 GByte of 
memory. As Fig. 1 shows, even these data are insufficient to extrapolate to the true asymptotic value. Numerical 
calculations that stop at even smaller values of n may easily misguide an extrapolation to c k (oo). 

In order to evaluate c k (n) for much larger values of n, we resort to another asymptotic result. In 1940, Mahler [3] 
showed that 

a k (n) = e**W V = e**< n > S k (n) (17) 
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Figure 2: Numerical evaluation of 



n 


C2(n) 


Sk(kn) Inn 
Sfc(n) n 


10 1 


1.49668 


1.50889 


10 2 


1.65470 


1.65496 


10 3 


1.65791 


1.65779 


10 4 


1.63881 


1.63876 


10 5 


1.61782 


1.61780 


10 6 


1.59883 


1.59882 


10 7 


1.58237 


1.58237 


10 8 


1.56822 


1.56822 


10 9 


1.55600 


1.55600 



Table I: Exact evaluation of C2(n) versus evaluation of the series. 



where </>fc(n) = 0{\). The idea is to replace the numerical evaluation of afc(n) by the numerical evaluation of the sum 
Sk(n). Note that S^{n) can be evaluated for very large values of n using a computer algebra system. A discrepancy 
in this approach arises from the unknown function <pk- Albeit asymptotically bounded, it can introduce large errors 
for finite values of n. 

It was already noticed by Froberg [5] , that <f>k oscillates with a small (and decaying) amplitude around a constant 
value. We used our extensive data for ak{n) to look more closely at 



: (n) = In a k (n) - ln^ 



(18) 



As can be seen from Figure 2, the amplitude of (f>2 is smaller than 10 4 for n > 10 , and it is slowly, but monotonically 
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l«j(2)| 


l«i(3)| 


1 


7.36616 ■ 10" 


7 


1.15010 • 10 -4 


2 


4.63909 ■ 10" 


13 


1.45894 ■ 10" 8 


3 


2.64857 ■ 10" 


19 


2.10798 • 10~ 12 


4 


1.29245 ■ 10" 


25 


1.45000 • 10~ 16 



Table II: Taylor coefficients (8) for 4>k and ipk- 
decaying. The constant around which (f>k oscillates will cancel in the ratio 

a k (n) S k (n) 

Hence the error is bounded by the small amplitude. This is confirmed by the numerical data, see Table I. Even for 
n = 100, the error in C2(n) is only in the fourth decimal. 
This observation tells us that we can use 

S k (kn) Inn 

b k(n) = -g-f-r: (19) 

Sk(n) n 

as an excellent approximation to Cfc(n). Since &&(n) can be evaluated for very large values of n, like n = 2 1000 and 
beyond, we can use bk(n) to bridge the gap between the numerically accessible Cfc(n) and Cfc(oo) (Figure 1). 



IV. ASYMPTOTICS RELOADED 



The results of dc Bruijn (6) and Mahler (17) have to match, i.e., we know that 4>k(n) + hi Sk(n) equals the right 
hand side of (6). A saddle-point expansion of Sk(n) (see (36) in the Appendix) reveals that the leading terms of 
hi Sk(n) equal the leading terms in (6). The remaining terms yield 

, . / / n \\ / lnlnfcN , 1 , , (lnfc + 21nlnA:) 2 /ln 2 lnn\ . . 



Jog fe n> 

In particular, we see that asymptotically <pk oscillates around a value 

lnlnfc\ , 1 . . (lnfc + 21nlnfc) 2 
1 + _ j In* - - hi(2n) - — + a (k) , (21) 

with ao(k) from (9). For k = 2, this constant is —0.079793025... (sec Figure 2), in perfect agreement with the 
numerical results of Froberg [5] . 

The asymptotic amplitude of 4>k is very small, as can be seen by evaluating the coefficients (8), see Table II. Hence 
we know that the oscillation in Figure 2 will eventually decay to an amplitude of size 10~ 6 . We have calculated a few 
more minima and maxima of 4>2 to check this decay. Figure 3 shows the result. The extrapolation of the numerical 
data gives very accurate result for the constant —0.079793025 ... as well as the right order of magnitude (10 -6 ) of 
the remanent amplitude. 



V. ACKNOWLEDGMENTS 



We thank Sebastian Mingramm for providing us with the numerical values for 0.2 (n) for n > 10 . SB thanks 
the Institut fur Thcorctische Physik at Otto-von-Guericke University in Magdeburg for its hospitality during the 
preparation of this manuscript and gratefully acknowledges support from a Fulbright-Kommision grant and from the 
U.S. National Science Foundation through grant number DMR-0812204. 



VI. APPENDIX 



Here, we evaluate the asymptotic behavior of the sum 



g fc(") = E— (»-°°) ( 22 ) 

i=or- k 2 
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Figure 3: Decay of the amplitude of <j>2- The scale of the abscissa is chosen to match 1 over the period of ipk- The numerical 
fits for the minima and the maxima are three parameter least square fits of fio + ^ix^ 2 . 



using a saddle-point expansion. Following Ref. [6] (pp. 304), we define <£>_,- = \naj for the summand a,j. The finite- 
difference condition = — = determines the maxima, i. e. we need to find the jo-te rm (s) of the sum 



with ctj /aj -i ~ 1. Applied to Eq. (22), we obtain 



77 [jo + Aj - 1 + log fc (jo + A,/)] ~ 1 



where we use the abbreviations 



V 



log fc n 



Ink 

Inn 



(23) 



(24) 



and a non- integer offset < Aj < 1 on the integer location jo of the saddlepoint, which we will need to attain the 
continuum limit for this n-dependent ("moving") saddle point [6]. For n — > oo, rj — > and in that limit we find from 
Eq. (23) for the saddle-point location by peeling off layer- by-layer: 



1 . , , , A 1 + lp gfc V , 2 1 + lQ gfc " 



inT + 1 + log ^ 



1 + logfc V 
61nfc 
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(25) 
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As jo 3> 1 f° r V ~~ > 0; we can expand <$>j for large arguments: 
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12j 



where we have used the Stirling expansion for the factorial to the necessary order. 

At the (unique) maximum of <frj we set j ~ j + t and expand here only to quadratic (Gaussian) order in t:[7] 
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(27) 



Note that the linear term in t only vanishes (indicating a symmetric maximum) after we insert the moving saddle-point 
in Eq. (25) and Aj is fixed: 



io+t 



ln(fc) 2 In (77) + In (fc) + 2 
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2?7 



In 2 fo) 

21n(fc) 



+ ln(r / )--ln(2 7 r) + 



3 ln(fc) 
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12 



lnry\ 
Ink) 



6- 



In 7/ 
lnfc 



lnfc 



+0 (r? 3 In 4 n) 



3 + 4^ + 2^ 
lnfc ln 2 fc 



7 



- (In (fc) + 77) (1 - 2A,) + ? ? 2 - " o A ^ + o A 



(In (fc) + r/) (1 - 2Aj) + 77 A 



12 2 J ' 2 1 
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t 2 + (1ft 3 ) . 



(28) 

(29) 
(30) 

(31) 



Terms of orders such as 7i 2 i 3 will not contribute at order rj 2 as they are at leading order asymmetric in t in the ensuing 
Gaussian integration. To that effect, we symmetrize the saddle point to order rj 1 with the choice of 

o2 



A,- = 



1 



+ 0(r/ 3 In 77), 



2 241nfc 

which also impacts constant or smaller terms in 77 in Eq. (30). The Gaussian integration then yields 
Sfc(n) 



e * 3o+t ' 



dt exp { -- 



(77 < e < 1) 
1 



(32) 



(33) 



ln(fc) + 77-r / 2 (l + I -^l 
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t 2 + 0(77 2 t 3 )|, 



r dt [l + 0(7 ? 2 t 3 )]exp{-i 
J —00 I 



1 + InT 



(34) 



1 



1 



exp <j $ jo + - In (2tt) - - In 



1 + InT 



+ O (t? 3 In 2 77) 



Note that the In (27r)-terms cancel. Hence, we finally obtain 
lnfc 2 In 77 + lnfc + 2 

v + 
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-v 
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(35) 
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or in terms of powers of Inn directly: 
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(36) 



In Fig. 4 we plot a sequence of approximants to the numerically exact evaluation of the sum in Eq. (22), which 
prove to approximate with an error of the indicated order. 
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(Concerned with sequence A033485 , A000123 , A005704 , A005705 , and A0057Q6.) 
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Figure 4: Plot of In Sk (n) / In 2 n vs 77 = lnfc/lnn for k — 2 (top set) and 3 (bottom set). The dotted horizontal lines specify 
the asymptotic limits, 2 \n.2 = 0.7213 . . . for k — 2 and 21 ^ 3 = 0.4551 . . . for k = 3. The thick black line is obtained from the 
numerically exact evaluation of Eq. (22), and the shaded, dashed lines correspond to the asymptotic expression in Eq. (35), 
evaluated to the indicated order in 77. 



